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We propose in this article an unambiguous definition of the local density of electromagnetic states 
Ohi (LDOS) in a vacuum near an interface in an equilibrium situation at temperature T. We show that 

the LDOS depends only on the electric field Green function of the system but does not reduce in 
general to the trace of its imaginary part as often used in the literature. We illustrate this result by 
^ — , a study of the LDOS variations with the distance to an interface and point out deviations from the 

standard definition. We show nevertheless that this definition remains correct at frequencies close 
to the material resonances such as surface polaritons. We also study the feasability of detecting 
such a LDOS with apetureless SNOM techniques. We first show that a thermal near-field emission 
spectrum above a sample should be detectable and that this measurement could give access to the 
electromagnetic LDOS. It is further shown that the apertureless SNOM is the optical analog of the 
scanning tunneling microscope which is known to detect the electronic LDOS. We also discuss some 
^ , recent SNOM experiments aimed at detecting the electromagnetic LDOS. 

PACS numbers: 03.50.De, 07.79. Fc, 44.40.+a, 73.20.Mf 
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The density of states (DOS) is a fundamental quantity from which many macroscopic quantities can be derived. 
_ 1 ■ Indeed, once the DOS is known, the partition function can be computed yielding the free energy of the system. It 
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I. INTRODUCTION 



follows that the heat capacity, forces, etc can be derived. A well-known example of a macroscopic quantity that 
follows immediately from the knowledge of the electromagnetic DOS is the Casimir forcei*^. Other examples are 
I ' shear forces 3 - and heat transfer— between two semi-infinite dielectrics. Recently, it has been shown that unexpected 
coherence properties of thermal emission at short distances from an interface separating vacuum from a polar material 
are due to the contribution to the density of states of resonant surface waves^. It has also been shown that the Casimir 
force can be interpreted as essentially due to the surface waves contribution to the DOSA*£. 

Calculating and measuring the local density of states (LDOS) in the vicinity of an interface separating a real 
material from a vacuum is therefore necessary to understand many problems currently studied. The density of states 
is usually derived from the Green function of the system by taking the imaginary part of the Green's function^. 
• i— i ' In solid-state physics, the electronic local density of states at the Fermi energy at the surface of a metal can be 
measured with a Scanning Tunneling Microscope (STM)£. This has been proved by several experiments, in particular 
(— | E the so-called quantum corral experiments^. Although one can formally generalize the definition of the electromagnetic 
O ' LDOS by using the trace of the imaginary part of the Green's tensor—, it turns out that this definition does not yield 
■ the correct equilibrium electromagnetic energy density. 

Recently, it has been shown theoretically that the STM and the Scanning Near-field Optical Microscope (SNOM) 
have strong analogies^. More precisely, in the weak tip-sample coupling limit, it was demonstrated that a unified 
formalism can be used to relate the STM signal to the electronic LDOS and the SNOM signal to electromagnetic 
LDOS. SNOM instruments^ have been used to perform different kinds of emission spectroscopy, such as lumines- 
cence^, Raman spectroscopj*i^ or two-photon fluorescence— For detection of infrared light, apertureless techniques^ 
have shown their reliability for imaging— as well as for vibrational spectroscopy on molecules^. Moreover, recent 
calculations and experiments have shown that an optical analog of the quantum corral could be designed, and that the 
measured SNOM images on such a structure present strong similarities with the calculated electromagnetic LDOS 22 23 . 
These results suggest that the electromagnetic LDOS could be directly measured with a SNOM. 

The purpose of this article is to show how the electromagnetic LDOS can be related to the electric Green-function, 
and to discuss possible measurements of the LDOS in SNOM. We first introduce a general definition of the electromag- 
netic LDOS in a vacuum in presence of materials, possibly lossy objects. Then, we show that under some well-defined 
circumstances, the LDOS is proportional to the imaginary part of the trace of the electrical Green function. The 
results are illustrated by calculating the LDOS above a metal surface. We show next that the signal detected with a 
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SNOM measuring the thermally emitted field near a heated body is closely related to the LDOS and conclude that 
the natural experiment to detect the LDOS is to perform a near-field thermal emission spectrum. We discuss the 
influence of the tip shape. We also discuss whether standard SNOM measurements using an external illumination can 
detect the electromagnetic LDOS^SiS 3 .. 

II. LOCAL DENSITY OF ELECTROMAGNETIC STATES IN A VACUUM 

As pointed out in the introduction, the LDOS is often defined as being the imaginary part of the trace of the 
electric-field Green dyadic. This approach seems to give a correct description in some cases^^ 3 -, but to our knowledge 
this definition has never been derived properly for electromagnetic fields in a general system that includes an arbitrary 
distribution of matter with possible losses. The aim of this section is to propose an unambiguous definition of the 
LDOS. 

Let us consider a system at equilibrium temperature T. Using statistical physics, we write the electromagnetic 
energy U(ui) at a given positive frequency to, as the product of the DOS by the mean energy of a state at temperature 
T, so that 

U ™ = ex P (W)-l (1) 

where 2irh is Planck's constant and ks is Boltzmann's constant. We can now introduce^ a local density of states by 
starting with the local density of electromagnetic energy energy U{v,u) at a given point r in space, and at a given 
frequency 10. This can be written by definition of the LDOS p(r,u) as 

U(r, lo) = p(r, w) (2) 
ex-p(nuj/ kb-1 ) — 1 

The density of electromagnetic energy is the sum of the electric energy and of the magnetic energy. At equilibrium, 
it can be calculated using the system Green's function and the fluctuation-dissipation theorem. Let us introduce the 
electric and magnetic field correlation functions for a stationnary system 

E l3 (v,v',t-t') = -L j d^e^vy ,u)e-^-^ = (E l {r,t)E*(v',t 1 )) (3) 

««(r,r , ,t-0 - ^ Jd U 7U i (Tj,u>)e-™W = (H i {T,t)H;{i>,1?)) (4) 

Note that here, the integration over u> goes from — oo to oo. If j(r) is the electric current density in the system, the 

<-,E 

electric field reads E(r, to) — ipou> J G (r, r', CL>).j(r')(i 3 r'. In the same way, the magnetic field is related to the density 

of magnetic currents m(r) by H(r, oj) = jG (r, r', cj)m(r')(i 3 r'. In these two expressions, G and G are the 
dyadic Green functions of the electric and magnetic field, respectively. The fluctuation-dissipation theorem yields^ 
that 



Tito eQiu 



= f „^..r^ ^ ImGf,{v,v'^) (5) 



[exp(huj/k B T) - 1] 2tt 

TlOJ 6qU. 

[exp(huj/k B T) - 1] ~2n 



If one considers only the positive frequencies U(r, uo) — 4 x e /2 J2i=i 3 £a{ r i r ! w ) + Mo/2 J2i=i 3 T~tu( r , r , w ) 
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[exp{Tiuj/kB-l ) — 1J 7r<r 



G (r,r,w) + G (r,r,w) 



so that 



(7) 



It is important to note that the magnetic-field Green function and the electric-Green function are not independent. 
In fact, one has 

2 H E 

-tG (r,r» = [V r x].G (r,r', «). [V r -x] (8) 
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A comparison of Eqs. J5J and Q shows that the LDOS of the electromagetic field reads 

- f(G E ) (9) 

in which /(G ) is an operator which will be discussed more precisely in the next section. 

III. DISCUSSION 

The goal of this section is to study the LDOS behavior for some well-characterised physical situations, based on 
the result in Eq. ©. 

A. Vacuum 

In a vacuum, the imaginary part of the trace of the electric and magnetic field Green functions are equal. Indeed, 
the electric and magnetic field Green functions obeys the same equations and have the same boundary conditions in 
this case (radiation condition at infinity). In a vacuum, the LDOS is thus obtained by considering the electric field 
contribution only, and multiplying the result by a factor of two. One recovers the familiar result 

uj 2 

Pv(r, w) = p v (uf) = -3-3 (10) 
which shows in particular that the LDOS is homogeneous and isotropic. 



p(r,uj) = — ^ImTr 

7T(r 



G (r,r,w) + G (r,r,w) 



B. Plane interface 



Let us consider a plane interface separating a vacuum (medium 1, corresponding to the upper half-space) from 
a semi- infinite material (medium 2, corresponding to the lower half-space) characterised by its complex dielectric 
constant 62 (u>) (the material is assumed to be linear, isotropic and non-magnetic). Inserting the expressions of the 
electric and magnetic- field Green functions for this geometry^S into equation @, one finds the expression of the LDOS 
at a given frequency and at a given height z above the interface in vacuum. In this situation, the magnetic and electric 
Green functions are not the same. This is due to the boundary counditions at the interface which are different for the 
electric and magnetic fields. In order to discuss the origin of the different contributions to the LDOS, we define and 
calculate an electric LDOS (p E (z, ui)) due to the electric-field Green function only, and a magnetic LDOS (p H (z, u>)) 
due to the magnetic-field Green function only. The total LDOS p E (z,uj) = p E (z,uj) + p H (z,oj) has a clear physical 
meaning unlike p E (z,uj) and p H (z,uj). Note that p E (z,u>) is the quantity which is usually calculated and considered 
to be the true LDOS. In the geometry considered here, the expression of the electric LDOS is 21 



P E (Z,LU) 



p v (uj) 



1 ndn 

n P 



2 + Re 



D 2ipujz/c 



Re 



(r? 2 



y 2ipujzfe 



(2k 2 - 1) 



[ ^[Im{rl 2 )+(2n 2 -l)l m {rl 2 )} 



-2\p\ujz I c 



(ii) 



This expression is actually a summation over all possible plane waves with wave number k = ui/c(k,p) where p = 
y/1 — k 2 if k < 1 and p = i\J n 2 — 1 if k > 1. rf 2 and r\ 2 are the Fresnel reflection factors between media 1 and 2 in 
s and p polarisations, respectively, for a parallel wave vector wk/o2&. < k < 1 corresponds to propagating waves 
whereas k > 1 corresponds to evanescent waves. A similar expression for the magnetic LDOS can be obtained : 

p v (ui) ( f 1 ndn r_ 



p H (z,u J ) = 



+ 



-1 UO 
ndn 



P 



Re 



(rf 2 e 2i ^ /c ) + Re (r^ 2 e 2 ^ /c ) (2k 2 - l) 



\P\ 



[Im(r p 12 ) + (2k 2 - 1) Im(r s l2 )] e - 2 H wz / c 



(12) 
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FIG. 1: LDOS versus frequency at different heights above a semi-infinite sample of aluminum. 



Adding the electric and magnetic contributions yields the total LDOS 

lo P 

p0 ° K 3 dn 



pM = P^tLL | J + (Re (r{ 2 e 2 ^' c ) + Re (r[ 2 e 2 ^/ c ))~ 

[I m {rl 2 )+Im{rl 2 )]e-^/ c \ (13) 



1 \p\ 

It is important to note that the electric and magnetic LDOS have similar expressions, but are in general not equal. 
The expression of p H (r) is obtained by exchanging the s and p polarisations in the expression of p E {v). As a result, 
the two polarisations have a symmetric role in the expression of the total LDOS p(r). 

The vacuum situation can be recovered from the previous expression by setting the values of the Fresnel reflection 
factors to zero. The same result is also obtained by taking the LDOS at large distance from the interface, i.e, for z ^ A 
where A = 2itc/lu is the wavelength. This means that at large distances, the interface does not perturb the density 
of electromagnetic states. In fact, e~ 2 l p l a,2 / c becomes negligible for the evanescent waves and e 2lpujz ^ c is a rapidly 
oscillating function for the propagating waves when integrating over k. The result is that all the terms containing 
exponential do not contribute to the integral giving the LDOS in the vacuum. 

Conversely, at short distance from the interface, p(r, u) is drastically modified compared to its free-space value. 
Equations lll|) - (|13fl show that the Fresnel coefficients and therefore the nature of the material play a crucial role 
in this modification. For example, as pointed out by Agarwal 2 ^, in the case of a perfectly conducting surface, the 
contribution of the electric and magnetic LDOS vanish, except for their free-space contribution. In this particular 
case, one also retrieves the vacuum result. 

We now focus our attention to real materials like metals and dielectrics. We first calculate p(cu) for aluminum at 
different heights. Aluminum is a metal whose dielectric constant is well described by a Drude model for near-UV, 
visible and near-IR frequencies 2 ^: 

e(w) = 1 - (14) 

with Lo p = 1.747 10 16 rad.s^ 1 and 7 = 7.596 10 13 rad.s" 1 . We plotted in FigJUthe LDOS p(r,u) in the near UV-near 
IR frequency domain at four different heights. We first note that the LDOS increases drastically when the distance 
to the material is reduced. As discussed in the previous paragraph, at large distance from the material, one retrieves 
the vacuum density of states. Note that at a given distance, it is always possible to find a sufficiently high frequency 
for which the corresponding wavelength is small compared to the distance so that a far-field situation is retrieved. 
When the distance to the material is reduced, additional modes are present: these are the evanescent modes that 
are confined close to the interface and that cannot be seen in the far field. Moreover, aluminum exhibits a resonance 
around uj = w p /v2. Below this frequency, the material supports resonant surface waves (surface-plasmon polaritons). 
Additional modes are therefore seen in the near field. This produces an increase of the LDOS close to the interface. 
The enhancement is particularly important at the resonant frequency which corresponds to Re[e(uj)] = — 1. This 
behavior is analogous to that previously described in Ref. 24 for a SiC surface supporting surface-phonon polaritons. 
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FIG. 2: Density of states contributions due to the propagating and evanescent waves compared to the toatal density of states 
and the vacuum density of states. These quantities are calculated above an aluminum sample at a distance of 10 nm. 




FIG. 3: LDOS at a distance z — Wnm above a semi-infinite aluminum sample. Comparison with p E {uj) and p H (w). 



Also note that in the low frequency regime, the LDOS increases. Finally, FigQ] shows that it is possible to have a 
LDOS smaller than that of vacuum at some particular distances and frequencies. Figure |21 shows the propagating and 
evanescent waves contributions to the LDOS above an aluminum sample at a distance of 10 nm. The propagating 
contribution is very similar to that of the vacuum LDOS. As expected, the evanescent contribution dominates at 
low- frequency and around the surface-plasmon polariton resonance), where pure near-field contributions dominates. 

We now turn to the comparison of p{z, u) with the usual definition often encountered in the literature, which 
corresponds to p E (z,oj). We plot in Fig[3|p, p E and p H above an aluminum surface at a distance z — lOnm. In this 
figure, it is possible to identify three different domains for the LDOS behaviour. We note again that in the far-field 
situation (corresponding here to high frequencies i.e. \/2ir <C z), the LDOS reduces to the vacuum situation. In 
this case p(z,oj) = 2p E (z,ui) — 2p H (z,ui). Around the resonance, the LDOS is dominated by the electric-field Green 
contribution. Conversely, at low frequencies, p H (z,oj) dominates. Thus, FigO shows that we have to be very careful 
when using the expression p(z) = p E (z,uj). Above aluminum and at a distance z — lOnm, this approximation is only 
valid on a small range between u> — 10 16 rads -1 and u> = 1.5 X 10 16 rarf.s _1 . 
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C. Asymptotic form of the LDOS in the near-field 

In order to get more physical insight, we have calculated the asymptotic LDOS behaviour in the three regimes 
mentioned above. As we have already seen, the far-field regime (X/2tt <C d) corresponds to the vacuum case. To 
study the near-field situation, we focus on the evanescent contribution as suggested by the results in Fig|3 When 
A = 2itc/uj 3> z, the exponential term exp(— \p\ujz/c) is small only for k ^> X/(Attz) 3> 1. In this (quasi-static) limit, 
the Fresnel reflection factors reduce to 



lim r 
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lim r^ 2 = 

K — >00 



4k 2 
e- 1 

e+1 



Asymptotically, the expressions of p E (z,u) and p (z, u>) are 

Pv 



p E (z,u) = 

l fc 

P H (z,uj) = p v 



l\ 2 Ak 3 z 3 



16k z A\e + l\ 2 k z 



(15) 
(16) 

(17) 
(18) 



At a distance z — lOnm above an aluminum surface, these asymptotic expressions matches almost perfectly with the 
full evanescent contributions (k > 1) of p E and p H . These expressions also show that for a given frequency, one can 
always find a distance to the interface z below which the dominant contribution to the LDOS will be the one due 
to the imaginary part of the electric-field Green function that varies like (koz)~ 3 . But for aluminum at a distance 
z = lOnm, this is not the case for all frequencies. As we mcntionncd before, this is only true around the resonance. 
For example for low frequencies, and for z — lOnm, the LDOS is actually dominated by p v e /(lGkoz). 



D. Spatial oscillations of the LDOS 



Let us now focus on the LDOS variations at a given frequency versus the distance z to the interface. There are 
essentially three regimes. First, for as discussed previously, at distances much larger than the wavelength the LDOS 
is given by the vacuum expression p v . The second regime is observed close to the interface where oscillations are 
observed. Indeed, at a given frequency, each incident plane wave on the interface can interfere with its reflected 
counterpart. This generate an interference pattern with a fringe spacing that depends on the angle and the frequency. 
Upon adding the contributions of all the plane waves over angles, the oscillating structure disappears except close 
to the interface. This leads to oscillations around distances on the order of the wavelength. This phenomenon is 
the electromagnetic analog of Friedel oscillations that can be observed in the electronic density of states near the 
interfaces'^. As soon as the distance becomes small compared to the wavelength, the phase factors exp(2ipu>z/c) in 
Eq. I|13|) are equal to unity. For a highly reflecting material, the real part of the reflecting coefficients are negative 
so that the LDOS decreases while approaching the surface. These two regimes are clearly observed for aluminum in 
Fig. El The third regime is observed at small distances as seen in Fig. 0] The evanescent contribution dominates and 
ultimately the LDOS always increases as 1/z 3 , following the behaviour found in Eq. <|17[) . This is the usual quasi-static 
contribution that is always found at short distance-^. At a frequency slightly smaller than the resonant frequency, 
surface waves are excited on the surface. These additional modes increase the LDOS according to an exponential law 
as seen in Fig. [S] a behavior which was already found for thermally emitted field o 5 i 26 . At low frequency, the LDOS 
dependance is given by Eq. l|18l) . The 1/z magnetic term dominates because the l/\e + 1| 2 takes large values. The 
1/z 3 contribution equals the 1/z contribution for distances much smaller than the nanometer scale, a distance for 
which the model is no longer valid. 

The main results of this section can be summarized as follows. The LDOS of the electromagnetic field can be 
unambiguously and properly defined from the local density of electromagnetic energy in a vacuum above a sample at 
temperature T in equilibrium. The LDOS can always be written as a function of the electric-field Green function only, 
but is in general not proportionnal to the trace of its imaginary part. An additional term proportional to the trace 
of the imaginary part of the magnetic-field Green function is present in the far-field and at low frequencies. At short 
distance from the surface of a material supporting surface modes (plasmon or phonon-polaritons), the LDOS presents 
a resonance at frequencies such that i2e[e(o;)] = — 1. Close to this resonance, the approximation p(z,u) = p E (z,u>) 
holds. In the next section, we discuss how the LDOS can be measured. 
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FIG. 4: LDOS versus the distance z from an aluminum-vacuum interface at the aluminum resonant frequency. 
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FIG. 5: LDOS versus the distance z from an aluminum-vacuum interface at frequency uj = 8 10 15 rad.s 1 . 



IV. MEASUREMENT OF THE LDOS 



A. Near-field thermal emission spectroscopy with an apertureless SNOM 



In this section we shall consider how the LDOS can be measured using a SNOM. We consider a frequency range 
where p is dominated by the electric contribution p E . We note that for an isotropic dipole, a lifetime measurement 
yields the LDOS as discussed by Wijnands et ali« 1Q . However, if the dipole has a fixed orientation x, the lifetime is 

proportionnal to G xx and not to the trace of G- In order to achieve a direct SNOM measurement of the LDOS, we 
have to fulfill two requirements. First, all the modes must be excited. The simplest way to achieve this is to use 
the thermally emitted radiation by a body at equilibrium. The second requirement is to have a detector with a flat 
response to all modes. To analyse this problem we use a formalism recently introduced. 

We consider a SNOM working in the detection mode, and detecting the electromagnetic field thermally emitted by 
a sample held at a temperature T. The system is depicted in Fig.El The microscope tip is scanned at close proximity 
of the interface separating the solid body from a vacuum. The signal is measured in the far field, by a point detector 
sensitive to the energy flux carried by the electromagnetic field. We assume that an analyzer is placed in front of the 
detector (polarized detection). The direction of polarization of the analyzer is along the direction of the vector j rec - 
If the solid angle dfl under which the detector is seen from the tip is small (a condition we assume for simplicity), the 
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FIG. 6: Scheme of a scanning near-field optical microscope measuring a thermally emitted field, (a): experimental situation, 
(b): reciprocal (fictitious) situation. 



signal (S(lu)) at the detector, at a given frequency w, reads 

(S{uj)) = ^ \E d (u)\ 2 r 2 dQ, (19) 

where eo is the permittivity of vacuum, c is the speed of light in vacuum, r is the distance between the tip and the 
detector, and Ed is the electric field at the position of the detector. Let us denote by E exp (experimental field) the 
thermal field, emitted by the sample, in the gap region between the tip and the sample. This field can be, in principle, 
calculated following the approach recently used in£*2i. For simplicity, we shall neglect the thermal emission from the 
tip itself (which is assumed to be cold) compared to that of the heated sample. But we do not need, at this stage, 
to assume a weak coupling between the tip and the sample. In particular, in the expressions derived in this section, 
the experimental field ~E exp is the field emitted by the sample alone, in the presence of the detecting tip. Following 
the approach o^, based on the reciprocity theorem of electromagnetism^, an exact relationship between the signal 
(S(lj)) and the experimental field ~Ei exp can be established. It can be shown that the signal is given by an overlapping 
integral. 

To proceed, one considers a fictitious situation in which the sample is removed, and a point source, represented 
by a monochromatic current j rec oscillating at frequency o>, is placed at the position of the detector (see Fig.HJb)). 
The orientation of this reciprocal source is chosen along the direction of polarization of the analyzer used in the 
experimental situation. The field created around the tip in this reciprocal situation is denoted by E rec . Using the 
reciprocity theorem, the field at the detector can be written^: 



E d M • irec = — I rec \j ' ; • E e:cp (R, z, u)dR (20) 



-2i r <9E ree (R, z,u) 
Houj J s dz 

where the integration is performed in a plane z = zq between the tip and the sample and R — (x,y) are the coordinates 
along this plane. Equation l|20[l connects the field above the surface E exp to the field in the detector Ed along the 
direction of the analyzer. Note that the reciprocal field E rec encodes all the information about the detection system 
(tip and collection optics). Reporting the expression of the field at the detector H2()J) in l|19(l . one finds the expression 
for the measured signal: 



(5( w )) = |2£ J J H ij (R,R',z,Lj)W ij (Ti,R',z,Lj)dIidR' 



(21) 



Equation 121|) establishes a linear relationship between the signal and the cross-spectral density tensor Wij of the 
electric field defined by 

(E e xp,i(R, z, Lo)E* exp j (H! , z, lo')) — Wij(Tl, R', z, lu)6(u> ~ lu') . (22) 

The response function Hij only depends on the detection system (in particular the tip geometry and composition), 
and is given by 

_ , (?E rec j(R, z, w) <3E rec j(TL',z, oj) ,„„. 

(R, R, u) = — q z — dz ' ^ ^ 
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The cross-spectral density tensor Wij describes the electric-field spatial correlation at a given frequency ui. For the 
thermal emission situation considered here, it depends only on the dielectric constant, on the geometry and on the 
temperature of the sample. 

Equation (|21[l is a general relationship between the signal and the cross-spectral density tensor. It is non-local and 
strongly polarization dependent. This shows that one do not measure in general a quantity which is proportional to 
Wkk{r 7 r, ui), and thus to p E (r,u)). Nevertheless Eq. (|2"T)l suggested that p E (r,uj) can be recovered if the response 
function Hij is localized. Indeed, in that case the signal is proportional to Wij (R, R, z, ui), thus to p E (v, ui). As shown 
in the next section, a dipole tip (small sphere) would exhibit such a response function. 



B. Detection of the LDOS by an ideal point-dipole probe 



Let us see what would be measured by an ideal probe consisting of a single electric dipole described by a polarizability 
a{ui). Note that such a probe was proposed as a model for the uncoated dielectric probe sometimes used in photon 
scanning tunneling microscopy (PSTM), and gives good qualitative prediction^. We assume that the thermally 
emitting medium occupies the half-space z < 0, and that the probe is placed at a point r t . As in the preceding 
section, the detector placed in the far field measures the field intensity at a given point Yd, through an analyser whose 
polarization direction is along the vector j rec . In this case, Eq. I|2UI) simplifies to read 

ui 1 expfi/clr^ — r t |) <-» . 

j rec ■ E d = a(uj) - — o ; ; j rec • h(u d ) ■ E exp (r t , w) (24) 

Aire 1 | r d- r t| 

where k = ui/c, Ud = (r^ — r t )/\rd — r t \ is the unit vector pointing from the probe towards the detector and 

h(ud) = I UdUd is the dyadic operator which projects a vector on the direction transverse to Ud, I being the unit 

dyadic operator. The dyadic h(ud) being symmetric, the scalar product in the right-hand side in Eq. I)24|) can be 

transformed using the equality j rec ■ h(ud) • E exp (r t , w) = E, exp (r t ,oj) ■ h(ud) • j rec . Finally, the signal at the detector 
writes 

(S) = \a(Lo)\ 2 ^dQj2A t A*W tJ (r u r t ,iu) (25) 

i,3 

where A = h(u<j) -j r ec is a vector depending only on the detection conditions (direction and polarization). Note that 
if jrec is transverse with respect to the direction u^, which is approximately the case in many experimental set-ups, 
then one simply has A = j rec . 

Equation (|25l) shows that with an ideal probe consisting of a signal dipole (with an isotropic polarizability a(uj), one 
locally measures the cross-spectral density tensor at the position r t of the tip. Nevertheless, polarization properties 
of the detection still exists so that the trace of Wij, and therefore p E (r, ui), is not directly measured. A possibility 
of measuring the trace would be to measure a signal (Si) in the direction normal to the surface with an unpolarized 
detection, and a signal (S2) in the direction parallel to the surface, with an analyzer in the vertical direction. (Si) would 
be a sum of the two signals obtained with j rec along the ir-direction and along the y-direction. (S2) would correspond 
to the signal measured with j rec along the z-direction. Using Eq. pSfl. we see that the signal (S) = (Si) + (S2) is 
proportional to the trace Wkk(^t,^t,uj), and thus to p E (r,uj). Measuring the thermal spectrum of emission with an 
apertureless SNOM which probe is dipolar is thus a natural way to achieve the measurement of p E (r,uj). Close to the 
material resonances, i.e in the frequency domain where p E (r,uj) matches p(r,u>), such a near-field thermal emission 
spectrum gives the electromagnetic LDOS. 



C. Analogy with scanning (electron) tunneling microscopy 

The result in this section shows that a SNOM measuring the thermally emitted field with a dipole probe (for 
example a sphere much smaller than the existing wavelengths) measures the electromagnetic LDOS of the sample in 
the frequency range situated around the resonant pulsation. As discussed above, the measured LDOS is that of the 
modes which can be excited in the thermal emission process in a cold vacuum. This result was obtained from Eq. 120|) 
assuming a weak tip-sample coupling, i.e., the experimental field is assumed to be the same with or without the tip. 

The same result could be obtained starting from the generalized Bardeen formula derived in ref^i. Using this 
formalism for a dipole probe, one also ends up with Eq. (|25p. which explicitly shows the linear relationship between 
the signal and p E (r,u>). This derivation is exactly the same as that used in the Tersoff and Haman theory of the 
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STM 8 . This theory showed, in the weak tip-sample coupling limit, that the electron-tunneling current measured in 
STM was proportional to the electronic LDOS of the sample, at the tip position, and at the Fermi energy. This 
result, although obtained under some approximations, was a breakthrough in understanding the STM signal. In the 
case of near-field optics, the present discussion, together with the use of the generalized Bardeen formulaii, shows 
that under similar approximations, a SNOM using an ideal dipole probe and measuring the field thermally emitted 
by the sample is the real optical analog to the electron STM. We believe that this situation provides for SNOM a 
great potential for local solid-surface spectroscopy, along the directions opened by STM. 

D. Could the LDOS be detected by standard SNOM techniques ? 

Before concluding, we will discuss the ability of standard SNOM techniques (by "standard" we mean techniques 
using laser-light illumination) to image the electromagnetic LDOS close to a sample. Recent experiments^ 3 - have 
shown that an illumination-mode SNOM using metal-coated tips and working in transmission produce images which 
reproduce calculated maps of p E (r,u>) (which is the adopted definition of the LDOS in this experimental work, see 
also—). We shall now show that this operating mode bear strong similarites to that corresponding to a SNOM 
working in collection mode, and measuring thermally emitted fields. This will explain why the images reproduce (at 
least qualitatively) the electric LDOS p E (r,ui). 

Let us first consider a collection-mode technique, in which the sample (assumed to be transparent) is illuminated 
in transmission by a monochromatic laser with frequency to, and the near-field light is collected by a local probe. If 
we assume the illuminating light to be spatially incoherent and isotropic in the lower half-space (with all incident 
directions included), then this illumination is similar to that produced by thermal fluctuations (except that only the 
modes corresponding to the frequency ui are actually excited). Note that this mode of illumination corresponds to 
that proposed in Ref. 33 , This similarity, together with the discussion in the precedding paragraph, allows to conclude 
that under these operating conditions, a collection- mode SNOM would produce images which closely resemble the 
electric LDOS p E {r,uj). 

We now turn to the discussion of images produced using an illumination-mode SNOM, as that used in RefiS 3 .. The 
use of the reciprocity theorem allows to derive an equivalence between illumination and collection-mode configurations, 
as shown in Ref^. Starting from the collection-mode instrument described above, the reciprocal illumination-mode 
configuration corresponds to a SNOM working in transmission, the light being collected by an integrating sphere over 
all possible transmission directions (including those below and above the critical angle). Under such conditions, the 
illumination-mode SNOM produces exactly the same image as the collection-mode SNOM using isotropic, spatially 
incoherent and monochromatic illumination. This explains why this instrument is able to produce images which 
closely follow the electric LDOS p E (r,u). Finally, note that in RefiS 3 ., the transmitted light is collected above the 
critical angle only, which in principle should be a drawback regarding the LDOS imaging. In these experiments, it 
seems that the interpretation of the images as maps of the electric LDOS remains nevertheless qualitatively correct, 
which shows that in this case, the main contribution to the LDOS comes from modes with wavevector corresponding 
to propagation directions above the critical angle. 

V. CONCLUSION 

In this paper, we have introduced a defintion of the electromagnetic LDOS p(r, ui). We have shown that it is fully 
determined by the electric-field Green function, but that in general it does not reduce to the trace of its imaginary 
part p E (r,co). We have studied the LDOS variations versus the distance to a material surface and have explicitly 
shown examples in which the LDOS deviates from p E (r,uj). Nevertheless, we have shown that around the material 
resonances (surface polaritons), the near-field LDOS reduces to p E (r,oj). Measuring the LDOS with an apertureless 
SNOM using a point-dipole tip should be feasible. The principle of the measurement is to record a near-field thermal 
emission spectrum. Under such condition, the instrument behaves as the optical analog of the STM, in the weak- 
coupling regime, which is known to measure the electronic LDOS on a metal surface. Finallly, we have discussed recent 
standard SNOM experiments in which the LDOS seems to be qualitatively measured. Using general arguments, we 
have discussed the relevance of such measurements and compared them to measurements based on thermal-emission 
spectroscopy. 
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APPENDIX A: CALCULATION OF THE FIELD AT THE DETECTOR FOR AN IDEAL POINT-DIPOLE 

PROBE 

Let e rec (K) and e exp (K.) be the tridimensional Fourier component of E rec (r) and E exp (r). In the configuration 
chosen in our problem the reciprocal field propagates to the negative z whereas the experimental fields propagates to 
the positive z. Thus 

E rec (r)= /e rec (K)exp[i(K.R- 7 (K)z)]dK (Al) 



E exp (r) = J e exp {K)exp[i(K.R + 1 (K)z)]dK (A2) 



where 7 (K) = ^Juj 2 /c 2 - K 2 . Putting Eqs. (|Alf> and l|A2jl into Eq.(J2DI) gives 



g 2 r 

E d (w) • jrec = / 7(K)e rec (-K) • e exp (K)dK 

^Mo J 



(A3) 



e exp (K) can be evaluated by calculating the field E exp (r). This last field is the field radiated by the reciprocal current 
j rec and diffused by the ideal probe. It can also be seen as the field radiated by the dipole induced at the position 
r t = (R t , z t ) of the probe. If p is the dipole induced at the position of the ideal probe, the reciprocal field at a position 
situated below the probe writes: 



E rec (r) = ^ w 2 G(r, r t ) • p 



2 f d2^gi[K.(R-R t )+ 7 (z t - 2 )] 



8tt 2 J 7 
where kg — uj 2 /c 2 . Comparing this expression and IjAlfl .thcn 



kk 

T 2 



(A4) 



erec(K) = ^^f ^^^^ ■ P (A5) 
where k~ = (K, —7). Furthermore, using the fact that h(k) = h(— k) and defining k + = (K,7) 

e rec (-K) = ^!e l ( K ' R '+^)E(k+) • P (A6) 

Let us denote E(jrec — > r t ) the field radiated by the reciprocal current in v t . The dipole induced then writes 
p = a(uj)e E(jrec — > r t ) and 

E(jrec^r t ) = -fl- h(u d ) ■ \ rec (A7) 

47T |r d -r t | 

Using the fact that for all dyadic A and for all vector a and b 

[A • a] • b = a • [A T • b] (A8) 

that h is a symetric dyadic (h = h ), that e exp {K) is transverse to the direction k + and the definition of E exp (r). 

. uj 2 exp(i/c|r,i — r t |) «-> 
jrec ■ E d = a(uj) - — j i i Jrec • h(u d ) ■ E, exp (r t ,Lu) (A9) 

AlTC Z \Vd - v t \ 
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